A NOTE ON THE MOTION OF SURFACES 



Robert I. McLachlan and Harvey Segur 

cn ■ 

C\ ; June 22, 1993 

as : 

(-H ' ABSTRACT. We study the motion of surfaces in an intrinsic formulation in which 

^ , the surface is described by its metric and curvature tensors. The evolution equa- 

5 ■ tions for the six quantities contained in these tensors are reduced in number in 

£f*^ ' two cases: (i) for arbitrary surfaces, we use principal coordinates to obtain two 

. equations for the two principal curvatures, highlighting the similarity with the 
equations of motion of a plane curve; and (ii) for surfaces with spatially constant 

t- H t negative curvature, we use parameterization by Tchebyshev nets to reduce to a 

^ i single evolution equation. We also obtain necessary and sufficient conditions for a 

CO ' surface to maintain spatially constant negative curvature as it moves. One choice 

^— 2 i for the surface's normal motion leads to the modified-Korteweg de Vries equation, 

■ the appearance of which is explained by connections to the AKNS hierarchy and 

' the motion of space curves. 

co ■ 

On i 
■ 

1. Introduction and equations of motion 

Studies of the motion of surfaces and interfaces occur in several scientific 
disciplines, including geometry [8,10,17], water waves [20], crystal growth [11,13], 
combustion [4], and gas dynamics [18]. Pelce [16] collected and reprinted an 
interesting collection of original articles in some of these areas. To simplify the 
analysis, these studies often consider only two-dimensional motion, thus reducing 
the problem to the motion of a curve [6], but this is only a special case of the 
usual situation. 

To study the evolution of two-dimensional surfaces arbitrarily embedded in 
three dimensions, a natural approach is to describe the surface purely intrinsi- 
cally, that is, in terms of the surface's metric and curvature tensors. In these 
variables, the equations of motion for the surface take the form of six coupled 
nonlinear evolution equations for the six independent components of the two 
tensors. These equations were first written down by Brower et al. [3], and later 
by Nakayama & Wadati [15]. 

Unfortunately, these equations are highly redundant. To see this, observe 
that if a surface that can be represented at a fixed time in the form z = z(x, y), 
then its motion can be described by a single evolution equation, for dz/dt. The 
purpose of this paper is to simplify the evolution equations obtained in [3] and 
[15], by reducing the number of unknowns. We obtain two simplifications, by 
choosing two different parameterizations of the surface. 

The first parameterization, in §2, is general: we parameterize a surface by its 
principal lines (i.e., lines tangent to the directions of principal curvature). This 
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reduces the system to four coupled equations, and, with a particular choice of 
parameterization along the principal lines, to two equations for the two curva- 
tures. If the surface happens to be a cylinder above a curve in a plane, and if 
the motion is two-dimensional, then one curvature vanishes and the equations 
reduce to the single equation obtained in [3,14] for the evolution of the curvature 
of a plane curve. 

The second parameterization, in §3, is more specialized: surfaces of constant 
negative curvature have an elegant parameterization by Tchebyshev nets. For 
surfaces that maintain constant negative curvature as they move, their motion is 
specified by a single evolution equation for the angle between the coordinate lines. 
As an example, we show that one choice for the normal velocity of the surface 
leads to the modified Korteweg-de Vries equation for the surface's evolution. 

For either set of coordinates discussed here, the key idea is to force the coor- 
dinate lines to maintain their defining property as the surface evolves, by adding 
an appropriate tangential velocity to the evolution. 

The first step is to derive the six coupled equations for the motion of a surface 
moving in R 3 . Local coordinates on the surface are (x\,X2) = x, time is t, and 
the position (in ]R 3 ) of a point on the surface is r(x, t). Later r (= dv/dt) will 
be specified, so one can think of x as a Lagrangian particle label on the surface. 
Define two tangent vectors and one normal vector: 

Ti X T 2 

r„ = r „ n = r , 



\Tl X T 2 | 

where greek indices range over n = 1, 2 with )M = tJ-. (Both the notation used 
here and the preliminary development follow Spivak [19, Chapter 2].) Note that 
n is a unit vector but the r M are not necessarily unit vectors. We have the metric 
g^y and curvature tensors (which define the first and second fundamental 
forms g^dx^dxy and h^dx^dx^, respectively): 

9^ = t^-t v , = n • t^ v = n • v^ v (1.1) 

As one moves along the surface (at a fixed time), the tangent and normal vectors 
change according to the Gauss- Weingarten equations, 

n,„ = -Txg^h^ 

where T x ^ u are the Christoffel symbols of the second kind defined by the met- 
ric 9fiu, and g^ v = {g^ v )~ l as usual. Repeated indices are summed on unless 
otherwise noted. 

The six quantities in (1.1) are not all independent. They are related by three 
consistency conditions for the PDE's (1.2) — the Gauss-Codazzi equations: 

#1212 = det/i (1.3a) 
Kx-^ = (1.36) 

where ;At denotes covariant derivative and -R1212 is the nontrivial element of 
the Riemann tensor. Note that (1.3b) gives independent information only for 
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(v, A, fi) = (1,1,2) or (2,2,1). Of the three remaining degrees of freedom, two 
are due to the arbitrary parameterization, leaving one — which we expect, since 
the surface could be written, e.g., r = (x, z(x),t). 

Therefore one might ask: What is to be gained by going to an "intrinsic" 
(i.e. g^ u , h^y) representation of the surface when eqs. (1.3) cannot usually 
be solved explicitly? Any surface motion can after all be formulated locally as 
i = /(x, t). Such an approach fails if the surface becomes vertical — an apparently 
artificial constraint — but the intrinsic formulation can also suffer from artificial 
coordinate singularities. One advantage is that in some applications, the motion 
of the surface normal to itself is given naturally in terms of the local curvatures 
[3,17] which are the eigenvalues of g^ x h\ v ; secondly, such an approach has been 
successful in the study of the motion of curves, as we now discuss. 

Hasimoto [9] showed that in the self-induction approximation, the motion of a 
3D vortex filament (which moves in the direction of its local binormal vector with 
speed equal to its local curvature) can be reduced to the nonlinear Schrodinger 
equation. Lamb [12] extended this to more general motions of curves, and found 
motions obeying the sine-Gordon and modified Korteweg-de Vries equations. All 
three of these equations are completely integrable [1] , and Goldstein and Petrich 
[7] showed how to specify infinitely many two-dimensional motions of curves 
in terms of integrable evolution equations. The occurrence of such integrable 
equations was shown in [14] to be due to the fact that the Serret-Frenet equations, 
(t = t/|t|, s =arclength) 



which are the analog of (1.2) for plane curves, are equivalent to the AKNS 



scattering problem [1] at zero eigenvalue. Thus evolutions I . I can be specified 



which give rise to integrable equations of the AKNS hierarchy. 

The equations for surfaces, (1.2), are more complicated than those for curves, 
(1.4), for two main reasons: (i) t in (1.33) is a unit vector, whereas the r M in 
(1.2) are not; and (ii) the coordinates, Xi , in (1.2) need not measure arclength. 
If we relax these two requirements in (1.4), then they become 



which are equivalent to (1.2) when \i = 1 only. 

There is another interesting connection between curves and surfaces which we 
note briefly. In [14] it is shown that the Serret-Frenet equations for space curves 
also reduce to the AKNS scattering problem at zero eigenvalue (i.e., they take 
the form (1.4) in appropriate variables). Considering a moving space curve as 
sweeping out a surface, we see that the equations analogous to Serret-Frenet for 
the surface, the Gauss- Weingarten equations (1.2), must also contain the AKNS 
scattering problem (in appropriate variables). This point of view explains the 
occurrence of integrable equations from the AKNS hierarchy in the equations 
defining particular surfaces, such as the sine-Gordon equation for surfaces of 
constant negative curvature. 




(1.4) 
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Motion of a surface can be prescribed in the form 

r = [/n + V% (1.5) 

where U and can be specified arbitrarily. This implies (differentiate r M = r ;At 
and n • r M = and use (1.2)): 

t, = + Vh^n + {-Uh\ + V\,)r x 

Now intrinsic evolution equations for can be derived in two ways. The 

first is to differentiate (1.1) with respect to t and to eliminate r M and n using 
(1.2, 1.3). This was first done by Brower [3] for = and by Nakayama et 
al. [15] for ^ 0. There is also an "intrinsic" derivation, namely requiring 
compatibility of the PDE's (1.2) and (1.6). This gives equations which must be 
solved for g^ u and h^, and the result is (of course) the same: 

9^ = -2hfj, v U + V mv + V v%lx 

V = U mv - h^h x „U + h^ x V\„ + h vX V\» + V x h^ x .„ 

Notice that does not affect the actual shape of the surface. The easiest 
way to see this is to note that the evolution in (1.5) is linear in U, V^, so we can 
consider U = 0; then the velocities just push particles around on a motionless 
surface. Notice also that if a coordinate system with some particular defining 
property is chosen initially, that property will be destroyed under general motion 
of the particles on the surface. Therefore we will use the tangential velocities V 1 
to "push" the coordinates back, so that they maintain their defining property as 
the surface moves. (In the next sections we show how to do this in some special 
cases.) Nakayama et al. [15] in addition allowed the coordinates to evolve on 
the surface, but this does not give any extra generality. 

2. Principal coordinates 

The principal curvatures (/ei, K2) and principal directions of a surface are the 
eigenvalues and eigenvectors of g^ x h Xv . If the coordinate lines are chosen to be 
tangent to the principal directions at each point on the surface, then g^ v and h^ v 
are diagonal, and we have principal coordinates. These can be defined globally 
on any surface, and are smooth except where K\ = K2 — umbilic points. (See [19, 
p. 288] for more information about umbilic points; they are generally isolated.) 
We define g\, g2, «i, and k 2 by 

_ ( e^ \ _ f e^Ki 

9lu/ ~ \ )' ^~\0 e*>K2 

Note that gi and do not form covariant vectors — they are just functions ap- 
pearing in the 2-tensors g^ and h^ v . 

In terms of these new variables, the Gauss-Codazzi equations (1.3b) now read 
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which we will use to simplify the subsequent equations. For completeness, we 
also give the Gauss- Weingarten formulae in these variables: 



















(=1 









=51-52 



-Kl 




(2.2) 



-k 2 

As the surface moves, the coordinate lines will not remain principal coordi- 
nates unless gi2 = /112 = 0. (This is where we begin to specify V 11 in order to 
maintain the desired coordinate system.) This leads to two equations in V^ 2 and 
V 2 -, whose solution is 

V] 2 = e~^f, V 2 1 = -e-^f (2.3) 

(U, 2 K 2 ),1 ~ (I7,l«l),2 



where 



/ 



(«2 - Kl) 2 



With these choices, (1.7) reduces to four coupled equations: 

9i = -2mU + V*g itj + 2V], (2.4a) 
= k 2 U + V'Kij + e-"/ 2 (e-"' 2 U,J,i + e-^'g ifV U tV /2 (2.46) 

= (1,2) or (2, 1); no sum on i, i'). These may be written as 
D gi 



Dt 



2V\ = -2 Ki U 



Dk- 

-^l = K 2 U + U, SiSi +e-^'g i , i ,U, i ,/2 



(2.5) 



where D/Dt is total derivative in the direction —V^dj (i.e. the right hand side 
of (2.5) gives the evolution at "unpushed" particle labels) and Sj is arclength 
along Xi, (i.e. s 1 = j Xl exp(±gi(x 1, x 2 )) dxi). 

As a special case of these equations, we note that if k 2 = d 2 = so that 
the surface only depends on xi, then / = in (2.3), we can take V 2 = 0, 
V 1 = W (xx) exp(-^gi) say, and we recover the equations for the motion of a 
plane curve [14]: 

g = 2W, S - 2kII k = k 2 U + U, ss + Wk, s 

Now we return to the general case, to reduce (2.4) to two coupled equations 
for the evolution of the two curvatures, «i and k 2 . These two equations, along 
with their compatibility condition, (1.3a), seem to be the simplest intrinsic for- 
mulation possible for a general surface. 

In (2.3), the equations for V^, there is still the freedom to specify two arbitrary 
functions, each of one variable. We can use this freedom to specify <?j, and hence 
the metric, completely. Here is one way to specify the metric: require 

g 1 (x 1 ,0,t)=g 2 {0,x 2 ,t) = V t, (2.6) 
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so that gn = g 22 = 1 along these two coordinate lines, which are both param- 
eterized by arclength (see Figure 1). Let t> M be two functions of integration in 
(2.3), i.e., 

V 1 (x 1 ,x 2 ) = + / exp(-g 1 (x 1 ,x 2 ))f{x 1 ,x 2 )dx2 + v 1 (x 1 ) 
Jo 

V 2 (x 1 ,x 2 ) = - ew(-92{xi,x 2 ))f(x 1 ,x 2 ) dx 1 + v 2 {x 2 ) 
Jo 



where we have suppressed the explicit t dependence. Then, from (2.4a) and 
choosing 1^(0) = 0, enforcing (2.6) requires 

^ = M-K^')| Ij(= o 
Now (2.1) can be integrated to define gi completely: 

9i= ■ — dx 2 , g 2 = - '- — dx 1 . 

Jo K 2 ~ Ki J K 2 - Ki 

Substituting these into (2.4b) yields two coupled equations for the evolution of 
the two curvatures, as desired. 

In the k 2 = d 2 = case, v[ = KiU and in this gauge, x\ is arclength, so we 
get W = f KiU ds as expected. 



3. Surfaces of constant negative curvature 

From the wealth of special surfaces and coordinate systems available, we con- 
sider a particular choice for which the equations of motion take a particularly 
simple form. At every point on a surface whose Gaussian (or total) curvature 
K is negative, the curvature tensor has two null vectors (i.e. vectors which 
annihilate h^, so v^v v h {lv = 0; see [19, p. 70]). An asymptotic line is a smooth 
curve on the surface whose tangent vector coincides with one of these null vec- 
tors at each point along the curve. If the Gaussian curvature is constant (in 
space) and negative, then there are two sets of asymptotic lines, and they form 
a Tchebyshev net, i.e., they can be parameterized by arclength everywhere [19, 
p. 368]. With these as coordinate lines, the metric and curvature tensors take 
the form 

(l cosuA / h 12 \ , . 

^= v cosc 1 J' h ^={h 12 J' (3 - 1} 

where u> measures the angle between the coordinate lines. Furthermore, the 
Gauss-Codazzi equations become: 

hi2 = V— K sin u, (3.2a) 

0^12 = — Ksmuj, (3.26) 

i.e., solutions of the sine-Gordon equation are in a local 1-1 correspondence with 
such surfaces. Unfortunately such coordinates may be only local [19, p. 373]. 
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The Gauss- Weingarten equations (1.2) also take the simple form 

(LU^i COt LU — Ci^i CSC (x> \ 

'_0 \_ y/^Ksmu T (3.3) 

\/—K cot LO —\/—K CSCLd J 

-(ln(/«i - k 2 )),i -(ln(/ci + k 2 )),i 

-k 1 k 2 /(k 1 - Ht 2 ) 

Ki + K 2 — («i - K 2 ) 

and similarly for T 2 , where 




and K\, k 2 are the principal curvatures V — -K"(cot cj =t cscw). 

The first step is to find those normal velocities U that keep the Gaussian 
curvature constant in space (but not necessarily in time). Substitute (3.1) into 
(1.7), and require both that the form of (3.1) be preserved (i.e., gn = g 22 = 
hn = h 22 = 0) and consistency (so that ui from (1.7a) equals that from (1.7b)). 
One finds in this way that U must satisfy a linear inhomogeneous equation: 

U, 12 = -KUCOSLU + (smuj)-^^K (3.4) 

and that the tangential components of velocity must satisfy: 

V\ v = -(KU smij+u^, ({-if U, v cos uj + U, i) + U, iH > smi^) / 2y/^K sm 2 uj (3.5) 

where i, i' = (1, 2) or (2, 1), and there is no sum on i'. There are no other con- 
straints. As mentioned earlier, only the normal component of velocity actually 
moves the surface, so we have established the following 

Proposition: Suppose that a surface whose Gaussian curvature K is constant in 
space and negative moves according to (1.5). It will maintain (spatially) constant 
negative curvature if and only if U satisfies (3.4). 

Note that if K = 0, then (3.4) is just the linearization of the sine-Gordon 
equation, (3.2b), so its translational symmetry immediately gives two solutions 
of (3.4): U = u t i and U = uj :2 . With the velocity components chosen to satisfy 
(3.4) and (3.5), the evolution of oj is given by 

+ ( v )2 + y 2 ,i) sin^ = 2^KU (3.6) 

where D/Dt was defined below (2.5). 

As an example, take K = and U = u t \ (=w ;Sl ). Eqs. (3.5) can be integrated 
to find 



V 1 = ((w,i) 2 + 2wncota;)/4^, V 2 = -w,n cscuj/2*J^K 

= I (6V=Ku tl + ((^ ;1 ) 3 + 2w,iu) /V^k) (3.7) 



so 

LO 
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which is the modified Korteweg-de Vries equation in the variable u \. The occur- 
rence of an equation from the AKNS hierarchy is not surprising, for the following 
reason. Let w(xi,0) be an arbitrary function. It can be integrated in x 2 under 
(3.2b) to sweep out a surface of constant negative curvature, or (because (3.7) 
does not depend explicitly on x 2 ) it can be integrated in t under (3.7). These 
two AKNS evolution equations commute — thus (3.7) obviously preserves the 
constant curvature property. In a sense, one should regard both x 2 and t as 
time-like variables. A similar instance of time-like evolution in spatial directions 
was reported in [2]. It is not known whether more complicated evolutions (e.g. 
U = + uj t2 ) can give rise to integrable equations. 

Now that it is known that this surface motion can be expressed as the mo- 
tion of a space curve it is fruitful to explore further. The coordinate lines are 
asymptotic and hence their principal unit normals lie in the surface (this can be 
seen from (3.3), because is proportional to the principal unit normal of the 
xi-curve). Furthermore for the particles in a coordinate line to trace out the 
coordinate grid properly it is easy to see (see Figure 2) that 

r ;2 = sin uj n c + cos ujt c = U c ii c + W c t c 

where the subscript c denotes properties of an xi -coordinate curve, k c is the 
curvature of the coordinate curve, and r c is its torsion. We know that this 
"motion" (i.e., motion in the ^-direction, as the curve sweeps out the initial 
surface) maintains the arc-length parameterization, which requires VF i = n c U =4> 
k c = Now take this curve motion and find n Cy2 using the space-curve 

evolution equations of [14]: the result is 

k C;2 = — t 2 sincj. 

Because u satisfies the sine-Gordon equation, we get r = y/—K, i.e., Tchebyshev 
coordinate lines in a surface of constant negative curvature have constant torsion 
[5] . The curve is evolving under the equation of motion 

r )2 = — sin / k c ds n c + cos / n c ds t c (3-8) 

We have in fact recovered the constant-torsion sine-Gordon curve equation of 
Lamb [12]. However, due to his inverse method (specifying t, n, b instead of r), 
he was unable to determine r for this evolution. In the forward direction it is 
easy to check that we have in fact the same evolution. 

(note: In the literature there are two distinct curve motions which obey the 
sine-Gordon equation. Let the curvature evolve by k c = QXJ . Lamb's evolution, 
given by (3.8), has QU = 0, and the curve would not move at all if r c = 0. 
Nakayama et al. [14] found a planar curve motion which obeys sine-Gordon: it 
has Q 2 U = 0.) 

As pointed out above, we can also think of an xi-coordinate line as moving 
in the t-direction, governed by the equation of motion (see Figure 2) 



r = Uh c + V 1 t c + F 2 (cosu;t c + sinu;n c ) 
= -K c b c + n cA n c /2V^K + K 2 c t c /4V^K 



(3.10) 
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The final cross-check is that a space curve with motion given by (3.10) has (from 
the curve equations in Nakayama et al. [14]) intrinsic evolution 

3 

k c = -V-Kk c , s + (3k^k C)S + 2k CjSSS )/4\/-K 

T C = 

agreeing with (3.7) when k c = — 
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Figure 1. Definition of the metric. 

An origin is fixed and two coordinate lines are chosen to be parameterized by 
arclength for all time; the metric at other points such as p is found by integrating 
the Gauss-Codazzi equations along a coordinate line. 



Figure 2. Tchebyshev net on a surface of constant negative curvature. 

All vectors shown are unit vectors, the normal n c to the coordinate curve lies in 

the surface, and the surface normal n is at right angles to it. 
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